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i^.n  experimental  investigation  was  cpnductet 
ter  mine  the  hot- gas  side  heat-transfer  rates  in  st  rocket  nozzle.  Transient  temperature 
measurements  were  made  at  five  axial  locations  in  a  copper  heat -sink  nozzle  having  ex¬ 
pansion  and  contraction  area  ratios  of  4.  64  with  gaseous  hydrogen  and  liquid  oxygen  as 
propellants.  The  nozzle  was  operated  over  a  range  of  chamber  pressures  from  150  to 
1000  pounds  per  square  inch  absolute  at  fuel  percentages  of  11  to  17,  with  most  measure¬ 
ments  made  at  approximately  15  percent  hydrogen  by  weight^' f 
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The  data  are  presented  in  a  correlation  equation  of  the  form  St  Pr  *  =  CRe^  *  , 


where  St  ,  Pr' 


and  Re*j  are  reference  Stanton,  Prandtl,  and  Reynolds  numbers,  re 


spectively,  and  C  is  a  constant.  Over  the  Reynolds  number  range  investigated,  the  re¬ 
sults  showed  that  no  one  constant  C  correlates  the  data  for  all  axial  locations;  how¬ 
ever,  the  data  did  correlate  well  when  each  station  was  handled  separately.  The  values 
of  C  varies  from  a  high  of  0.  026  in  the  chamber  to  a  low  of  0.  015  at  the  throat.  This 
throat  value  is  42  percent  lower  than  the  widely  used  value  of  C  =  0.  026. 
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INTRODUCTION 

(•  In  the  design  of  convectively  cooled  nozzles  for  high-performance  chemical  and  nu¬ 
clear  rocket  systems  the  high  heat  fluxes  and  associated  high  wall  temperatures  require 
accurate  knowledge  of  the  heat  flow  across  the  gas  boundary.  Present  design  procedures 
involve  gas-side  heat-transfer  correlations  which  have  not  been  well  substantiated  by 
rocket  heat-transfer  data.]  Rocket-nozzle  heat-transfer  investigations  are  reported  in 
references  1  to  5.  Although  the  data  of  references  1  to  5  appear  to  correlate  moderately 
well,  the  chamber  pressure  range,  propellant  combination,  nozzle  geometry,  and  in¬ 
strumentation  technique  limit  the  usefulness  of  the  data. 


I  It  was!  therefore^  deemed  desirable  to  conduct  a  gas-side  heat-transfer  experiment 
over  a  wide  range  of  chamber  pressures  with  emphasis  on  the  accurate  determination  of 
gas-side  heat-transfer  rates  in  the  regime  of  high  heat  flux.  Although  this  program  is 
directed  toward  the  solution  of  problems  associated  with  both  the  cooling  side  and  the 
hot-gas  side,  the  investigation  reported  herein  had  the  objective  of  providing  useful  hot- 
gas-side  design  information.  \  ()  v  <a 

The  rocket  nozzle  used'for  this-mvestigation  had  a  5 -inch- diameter  throat  with  a 
4.  64  expansion  and  contraction  area  ratio.  The  propellants  used  were  gaseous  hydrogen 
and  liquid  oxygen.  The  nozzle  wall  was  constructed  of  copper,  which  provided  heat-sink 
cooling  for  the  transiently  conducted  tests.  Local  heat-transfer  rates  were  determined 
at  five  axial  stations  through  the  nozzle  by  transient  temperature  measurements  using 
instrumented  copper  rods.  The  range  of  chamber  pressures  covered  was  approximately 
150  to  1000  pounds  per  square  inch  absolute.  The  range  of  mixture  ratios  investigated 
corresponded  to  11  to  17  percent  hydrogen  by  weight,  with  the  bulk  of  the  runs  being 
made  at  15  percent  hydrogen. 


APPARATUS  AND  EXPERIMENTAL  PROCEDURE 
Rocket  Design 


A  copper  heat-sink,  solid-wall,  hydrogen- oxygen  rocket  was  used  in  the  experiments 
to  obtain  transient  temperature  data  (fig.  1).  The  nozzle  had  a  4.  64  contraction  and  ex- 


Figure  1.  -  Copper  heat-sink  rocket  nozzle  with  coaxial  injector.  Dimensions  given  in  inches. 
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(a)  Injector  without  faceplate. 


(b)  Injector  with  copper  faceplate. 
Figure  2.  -  Coaxial  234  element  injector. 
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pansion  area  ratio.  The  length  of  the  combustion  chamber  from  the  injector  to  the  throat 
was  14.  5  inches.  The  characteristic  length  L*  was  54  inches.  The  throat  diameter 
was  5  inches.  The  rocket  developed  26  000  pounds  of  thrust  at  a  chamber  pressure  of 
900  pounds  per  square  inch  absolute. 

A  coaxial  injector  (figs.  1  and  2)  was  used  with  liquid  oxygen  and  gaseous  hydrogen 
as  the  propellants.  The  injector  had  234  injector  elements  uniformly  spaced  in  a  dish¬ 
shaped  faceplate  (2.  6  elements  per  square  inch).  Both  chamber  and  injector  were  de¬ 
signed  with  the  idea  of  obtaining  stable  engine  operation  and  as  uniform  a  distribution  of 
mass  flux  and  temperature  as  possible  to  eliminate  striation  effects. 


Instrumentation 


A  small  copper  plug  or  rod  was  devised  that,  when  inserted  in  the  wall  of  the  un- 
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cooled  rocket  engine,  approximated  the  model  of  a  one-dimensional  semi-infinite  slab. 

As  shown  in  figure  3,  holes  were  machined  normal  to  the  inside  wall  of  the  copper  nozzle 
at  each  of  the  seven  measuring  stations  (five  axial  locations  and  two  additional  circum¬ 
ferential  locations  at  the  throat).  Rods  of  pure,  oxygen-free  copper  with  diameters  of 
0.  25  inch  were  used.  The  rods  were  threaded  at  the  end  for  insertion  in  the  machined 
holes.  The  threaded  portion  was  machined  so  that  its  cross-sectional  area  closely  ap¬ 
proximated  that  of  the  0.  25-inch-diameter  rod.  This  threaded  portion  was  oxidized  be¬ 
fore  insertion  to  minimize  conduction  between  the  rod  and  the  chamber.  The  length  of 
the  rods  was  2.  5  inches,  and  the  thickness  of  the  chamber  wall  2  inches.  This  length 
was  picked  to  minimize  the  temperature  differential  between  the  heat-transfer  surface  of 
the  rod  and  the  surrounding  chamber  wall  because  of  the  larger  volume  of  a  radial  ele¬ 
ment  compared  to  the  volume  of  the  rod. 

Figure  3  also  shows  the  thermocouple  stations  on  the  copper  rod.  Six  Chromel- 
Alumel  thermocouples  were  used.  The  Chromel  and  Alumel  wires,  0.  008  inch  in  diam¬ 
eter,  were  peened  to  the  surface.  Figure  3  shows  a  detail  of  the  sealing  method  used  at 
the  end  of  the  rod.  The  cavity  surrounding  the  rod  was  pressurized  with  room- 
temperature  nitrogen  to  a  pressure  5  to  10  percent  below  that  of  the  local  static  pressure 
at  each  measuring  station  to  prevent  hot  gas  leakage  past  the  threaded  portion  of  the  rod. 
The  pressures  were  monitored  during  each  run  to  ensure  that  no  film  cooling  was  ob¬ 
tained  from  the  pressurizing  gas.  Cavity  pressure  remained  constant  with  time  after  the 
start  of  a  run. 


Data  Recording 

Propellant  flow,  chamber  pressure,  and  rod  temperatures  were  recorded  in  a  digi¬ 
tized  form  on  a  magnetic  tape  for  direct  entry  into  a  digital  computer.  The  digital  re¬ 
cording  system  had  basic  sampling  rates  of  4000  words  per  second.  Figure  4  (p.  6) 
shows  the  basic  block  of  data  as  taken  by  the  digitizer.  Most  data  parameters  were  ar¬ 
ranged  in  the  block  so  that,  with  a  minimum  of  five  blocks  of  data,  60- cycle  noise  could 
be  defined.  Then  a  smoothing  technique  was  used  to  eliminate  the  60- cycle  noise  and 
greatly  diminish  any  random  noise.  The  calculating  procedure  fit  a  smooth  curve  through 
15  readings  of  the  data  parameter,  and  then  one  report,  to  be  used  in  the  terminal  calcu¬ 
lations,  was  made  for  all  parameters  at  a  common  time.  This  reduced  the  amount  of 
terminal  calculations.  Chamber  pressure,  which  was  sampled  nine  times  in  each  block, 
was  used  as  a  triggering  device  for  starting  and  stopping  the  calculations. 
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Figure  4.  -  Digitizer  data  block. 


Test  Procedure 


The  rocket  engine  (fig.  5)  was  installed  on  a  test  stand  located  at  the  Lewis  Plum 
Brook  FacilS  Propellant  valves  for  controlling  gaseous  hydrogen  and  liquid  oxygen 
were  prepositioned  hire  the  run  and  opened  to  these  fixed  positions  durmg  the  run  to 
provide  the  desired  values  for  chamber  pressure  and  percent  of  hydrogem  Chamber 
pressures  were  varied  from  150  to  1000  pounds  per  square 

cent  hydrogen.  The  timing  and  sequencing  of  the  valves  were  adjusted  so  that  full 
her  pressure  was  obtained  in  0.  02  to  0.  06  second  for  a  typical  run  as  shown  m  figu 
This  step  function  in  chamber  pressure  or  driving  temperature  allowed  a  sxmpler  mat 
matical  model  to  be  used  for  obtaining  heat-transfer  coefficients. 


CALCULATION  PROCEDURE 
Constant  h  Method 

The  differential  equation  for  one-dimensional  heat  conduction  in  a  homogeneous 
medium  is 
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Chamber  pressure, 
Ib/sq  in.  abs 


Figure  5.  -  Copper  heat-sink  nozzle  installed  in  test  facility. 


Time,  sec 


Figure  6.  -  Typical  oscillograph  trace  of  chamber  pressure  against  time. 
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(1) 
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3X2 


pc 


ST 

St 


with  the  boundary  conditions 


-k^'h<Taw-T>  atx-° 


(2) 


k  —  =  0  at  X  =  if  =  °°  (3) 

3X 


The  initial  condition  is  T  =  Tq  at  t  =  0.  (Symbols  are  defined  in  the  appendix. )  The 
assumption  that  the  heat-transfer  rates  are  proportional  to  the  difference  between  the 
temperature  at  the  wall  and  the  temperature  of  the  surrounding  medium  (eq.  (2))  seems 
to  be  justified  because  the  heat  transfer  occurs  primarily  by  convection,  rather  than 
conduction  or  radiation.  Radiant  heat  transfer  is  estimated  to  be  an  order  of  magnitude 
less  than  that  due  to  convection  for  this  experiment.  It  is  therefore  neglected  in  the 
theoretical  treatment.  The  measurements,  however,  include  both  effects.  The  solution 
to  equation  (1)  is  given  in  closed  form  in  reference  6. 

Reference  7  has  put  this  solution  in  the  form 


^(Nx,7*)  =  erfc  — — - e 

2yfr* 

where 


erfc  z  =  1  — — 


while  the  dimensionless  variables  are 
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<p  = 


T  -  T, 


T  -  T 
aw 


0 


N, 


hX 

k 


T*  =  h^t 
kpc 

Equation  (4)  has  the  factor  2  added,  which  was  erroneously  omitted  in  reference  7. 

Equation  (4)  holds  until  such  time  as  the  temperature  starts  to  rise  at  the  position 
X  =  &  (the  thickness  of  the  slab).  For  shorter  times  the  temperature  distribution  given 
by  equation  (4)  is  approximately  represented  by  that  in  a  semi-infinite  solid,  the  accuracy 
increasing  as  X  =  £3?  approaches  the  hot  boundary  (£  -  0). 

When  the  temperature  starts  to  rise  at  X  =  if,  equation  (4)  fails  because  of  boundary 
condition  (eq.  (3)).  Since  longer  running  times  are  needed,  this  failure  of  the  solution 
can  be  remedied  by  allowing,  in  addition  to  the  heat  source  of  strength  h(T  -  T)  at 
X  =  0,  an  additional  heat  source  of  this  strength  at  X  =  2  if.  Reference  7  shows  that, 
with  two  identical  sources  located  symmetrically  about  X  =  <£,  the  conduction  of  heat  in 
the  region  0  <  X  <  if  is  governed  by  the  same  differential  equation  and  initial  boundary 
conditions  as  in  the  slab,  but  the  solution  given  applies  during  the  time  when  the  temper¬ 
ature  at  X  =  0  is  not  influenced  by  the  image  source  at  X  =  2  if.  The  temperature  of 
the  slab  is  therefore  accurately  represented  within  this  extended  time  period  by  the 
superimposed  effect  of  the  two  sources: 

^(NX5  T  )  =  ^(NX=|if ’  T*)  +  ^[NX=(2-£)i? ’ 

Equation  (5)  was  programed  on  a  digital  computer  so  that  an  iteration  process  found 
the  h  that  satisfied  the  measured  (p  for  a  given  location  |  and  time  r*.  The  last 
term,  <P[Nx=(2-£).2”  T*]>  was  a^so  reP°rted  so  that  cutoff  of  data  could  be  made  when  the 
temperature  rise  at  X  =  0  was  not  yet  influenced  by  the  image  source  at  X  =  2if.  Since 
the  engine  runs  were  all  necessarily  of  short  duration,  equation  (5)  was  sufficient  for 
reduction  of  the  data. 

The  properties  of  materials  (k,  p,  and  c)  may  vary  appreciably  with  temperature, 
and  a  problem  arises  concerning  which  temperature  in  the  material  should  be  used  in  the 
evaluation.  Reference  8  shows  that  the  material  properties  for  simple  metals  for  this 
type  of  computation  should  be  evaluated  at  a  temperature  that  is  about  one-fourth  of  the 
maximum  if  the  initial  temperature  is  zero,  or  at  (T^_0  -  TQ)/4  +  TQ  if  the  initial  tem- 
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Figure  7.  -  Proportionality  constant  and  nondimensional 
temperature  as  functions  of  time.  Distance  normal  to 
heated  surface,  0. 125  inch;  k/pc,  0. 001151  square  foot 
per  second. 


perature  is  nonzero.  The  tempera¬ 
ture  of  the  rod  at  £  =  0  was  approxi¬ 
mated  with  finite  differences  by  the 

equation 


TS=0=3TH>.05'2TH)-  l  +  3  TH>.2 

where  T^=0  Q5,  T^=0  v  and  T^=0  2 
are  measured  values. 

Figure  7  shows  the  results  of  a 
study  that  relates  the  error  in  a  non- 
dimensional  temperature  cp  to  the 
percent  error  in  the  nondimensional 
heat  transfer  hX/k.  The  nondimen¬ 
sional  temperature  cp  is  plotted 
against  time  for  various  hX/k  val¬ 
ues.  The  proportionality  constant 
G  from  the  equation 

hX/k 

is  plotted  against  time  for  the  same 
hX/k  values.  Since  the  limit  of 
running  time  for  a  heat- sink  engine 
is  a  constant  value  of  cp  (Tw  close 
to  melting  at  the  surface),  it  can  be 


seen  from  figure  7  that  G  remains  essentially  constant  with  hX/k.  Thus,  even  though 
the  length  of  running  time  goes  down  as  chamber  pressure  is  increased  (hX/k  is  in¬ 
creased),  the  accuracy  of  obtaining  h  remains  the  same  for  a  given  measurement  in¬ 
accuracy  dtp  and  a  constant  cp.  If  the  accuracy  in  obtaining  cp  is  estimated  to  be 
±5  percent  in  the  range  of  cp  =  0.  15,  the  percent  error  in  hX/k  is  ±8  (0.  15)  (0.  05)  (100) 


or  ±6.  0  percent,  as  shown  by  arrow  in  figure  7. 


Integration  Method 


An  integration  method  was 


also  used  to  check  the  results  from  the  constant  h  method 


described.  If  the  slab  is  considered  to  be  insulated  on  all  surfaces  except  the  heat- 
transfer  surface,  the  time  rate  of  change  of  the  integrated  heat  content  must  be  equal  to 
the  rate  of  heat  flow  through  the  surfaces: 


11  ^  AVQV  dX  =  !  —  (pAcT)x  dX  (6) 

A3tjJ  X  X  A  3t  ^ 


In  equation  (6)  the  heat  content  Qx  is  above  some  arbitrary  reference  level.  If  the  slab 
cross-sectional  area  is  constant  over  its  length,  the  area  term  does  not  appear  in  equa¬ 
tion  (6).  If  a  short  time  increment  after  the  start  of  steady-state  burning  is  considered, 
the  time  rate  of  change  of  integrated  heat  content  at  the  arithmetic  average  of  these  two 
times  is  equal  to  the  difference  in  heat  content  divided  by  the  time  interval.  The  heat- 
transfer  coefficient  can  be  obtained  from  the  following  equation: 


h0.  5(^+4)  "  7^ 

1  z  ^aw 


T  ) 
w' 


0.  5(t  i+tg) 


In  this  equation  the  pc  product  is  evaluated  at  the  measured  temperature  at  X. 


Combustion  Temperature 

Before  h  can  be  computed  by  the  methods  described  previously,  the  driving  tem¬ 
perature  (T  )  must  be  determined  throughout  the  nozzle;  to  determine  Taw  the  com¬ 
bustion  temperature  must  be  known.  Combustion  temperature  is  a  function  of  combustion 
chamber  pressure,  percent  fuel,  and  combustion  efficiency. 

Numerous  programs  are  presently  in  existence  for  the  calculation  of  equilibrium 
compositions  and  other  thermodynamic  properties  of  complex  chemical  systems.  The 
programs  of  reference  9  for  thermodynamic  properties  and  of  reference  10  for  transport 
properties  were  modified  and  simplified  for  a  gaseous-hydrogen  -  liquid-oxygen  system 
by  Frederic  N.  Goldberg  of  the  NASA  Lewis  Research  Center.  The  program  was  further 
modified  to  account  for  combustion  efficiency.  The  procedure  used  to  account  for  com¬ 
bustion  efficiency  and  arrive  at  the  proper  combustion  temperature  is  described  in  the 
paragraphs  that  follow. 

Under  the  conditions  of  one-dimensional  isentropic  equilibrium  flow,  the  specifica¬ 
tion  of  chamber  enthalpy  Hc  (for  a  given  mixture  ratio),  combustion  chamber  pressure 
p  and  throat  area  per  unit  weight  flow  (A/ w)^  may  lead  to  an  incompatability.  If  the 
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calculated  value  of  (A/w)mjn  corresponding  to  the  minimum  value  of  the  function 
A/w  =  f(P,  S  )  (where  S  is  the  entropy  of  the  chamber  gases  determined  by  H  and 
Pc)  differs  from  (A/w)^,  then  one  of  the  input  specifications  must  be  relaxed. 

Since  the  measurements  considered  to  be  most  reliable  were  those  of  combustion 
pressure  and  geometric  throat  area,  these  were  considered  to  be  invariant.  When  the 
calculated  value  (A/w)min  was  greater  than  (A/w)th,  an  iteration  was  used  to  determine 
a  chamber  enthalpy,  or  equivalently  a  combustion  temperature  Tc,  compatible  with  the 
measured  weight  flow  and  the  requirement  that  choking  occurs  at  the  geometric  throat 
according  to  the  following  algorithm: 


=  f(p,sK) 


tk+i  tk 

c  c 


sf+1  =  s(t*+1,p) 


where  the  K  =  0  conditions  correspond  to  those  of  specified  input  enthalpy.  Calculations 
were  terminated  when  there  was  sufficiently  close  agreement  between  the  calculated  and 
measured  throat  area  per  unit  weight  flow.  The  square  root  of  the  ratio  of  the  final  com- 
bustion  temperature  to  T£  was  used  as  a  measure  of  combustion  efficiency. 

If  the  initial  calculated  (A/w)min  was  less  than  that  measured,  the  enthalpy  was  left 
invariant.  (Note  that  in  this  case  the  preceding  method  would  indicate  an  efficiency 
greater  than  100  percent. )  A  new  weight  flow  was  then  calculated  with  the  combustion 
efficiency  assumed  to  be  100  percent: 


measured 


The  average  combustion  efficiency  obtained  from  all  the  test  data  was  98  percent  of 
theoretical  equilibrium. 


Transport  Properties 


Eckert  and  Drake,  in  a  summary  discussing  turbulent  boundary  layers  (ref.  11), 
recommend  that  the  transport  properties  appearing  in  the  correlation  equations  be  in¬ 
troduced  at  a  reference  enthalpy  given  by  the  equation 

H*  =  Hg  +  0.  5(HW  -  Hg)  +  0.  22(Pr*)1/3(Htot  -  Hg) 

instead  of  by  temperature  when  the  range  of  temperature  is  such  that  the  specific  heat 
varies  considerably  within  the  boundary  layer.  Therefore,  the  nondimensional  heat- 
transfer  parameters  were  computed  by  introducing  the  transport  properties  as  a  function 
of  H*  and  P  for  the  hydrogen- oxygen  system  described. 
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Figure  8.  -  Typical  profiles  of  wall  temperature,  heat  flow 
rate,  and  heat-transfer  coefficient  plotted  against  area 
ratio.  Chamber  pressure,  967  pounds  per  square  inch 
absolute;  fuel  mixture,  16. 5  percent  hydrogen;  time, 
0.462  second. 


RESULTS  AND  DISCUSSION 

Figure  8  is  a  plot  of  the  wall 
temperature,  heat  flux,  and  heat- 
transfer  coefficient  for  a  typical  run, 
for  which  the  chamber  pressure  was 
967  pounds  per  square  inch  absolute, 
the  hydrogen  concentration  was 
16.  5  percent,  and  the  time  0.  462 
second.  Figure  8  shows  that  the  heat 
flux  was  indeed  in  the  range  of  anti¬ 
cipated  nuclear-rocket  nozzles,  that 
is,  q  =  20  Btu  per  square  inch  per 

second,  and  T  =  2000°  R. 
w 

The  results  presented  in  this 
report  are  all  from  the  constant  h 
solution.  As  running  progressed, 
thermocouples  at  X/if  =  0.  6  and 
X/if  =  1.  0  were  lost  and  the  inte¬ 
gration  method  could  not  be  applied. 
For  some  of  the  earlier  tests,  how¬ 
ever,  results  from  the  integration 
method  agreed  with  the  constant  h 
solution  and  gave  additional  support 
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Figure  9.  -  Trend  comparison  of  theoretical  and  measured  temperature  transients  in  semi¬ 
infinite  solid.  (x  =  0.125  in.) 


to  this  method.  For  a  further  test  of  the  apparatus  and  technique  a  second  nozzle  fabri¬ 
cated  and  instrumented  in  the  same  manner  was  run  at  two  of  the  chamber  conditions. 
Data  reproducibility  was  shown  to  be  good;  that  is,  it  fell  within  the  spread  of  the  data 
and  thus  indicated  nothing  unique  about  one  particular  apparatus. 

Before  the  results  of  any  experiment  can  be  used,  the  validity  of  the  results  must  be 
established.  Figure  9  can  be  used  to  consider  the  results.  The  solid  lines  are  the  theo¬ 
retical  lines  from  equation  (4)  rewritten  as  follows  to  include  the  pertinent  variables 


for  hX/k  =  0.  05,  0.  075,  0.  1,  and  0.  125  and  X  =  0.  125  inch.  Experimental  data  from 
one  thermocouple  for  the  throat  station  are  shown  for  three  chamber  pressures.  The 
experimental  points  follow  the  contours  of  the  theoretical  lines  fairly  well;  thus  the  as¬ 
sumption  of  a  constant  h  solution  and  a  one-dimensional  model  seems  justified. 

The  data  for  the  range  of  chamber  pressures  from  150  to  1000  pounds  per  square 
inch  absolute  are  presented  in  figure  10  for  each  measuring  station  of  the  rocket  nozzle. 
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Stanton  number  X  (Prandtl  number)0- 1 ,  St  Pr 


Figure 
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Constant, 

C 
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deviation, 

<7 

-  x  100, 

c 

percent 
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c 

per  cent 

10(a) 

1 

0.  0257 

0.  0028 

10.9 
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2 

.0240 

.0039 

16.3 

.0077 

32. 1 

10(c) 

3 

.0148 

.  0020 

13.  5 

.0039 

26.3 

10(d) 

3 

.0151 

.0020 

13.2 

.0039 

25.8 

10(e) 

4 

.0153 

.0016 

10.4 

.0031 

20.3 

10(f) 

5 

.0188 

.0015 

7.  9 

.0030 

16.0 

(a)  Correlation  equation,  St  Pr  1 


4  3  2  1  2  3 

Area  ratio,  A/Av 

(b)  Correlation  equation,  St  Pr  '0- 7  =  CRejjb. 


Station  3  represents  the  average  values 

for  the  three  measuring  stations  at  the 

throat  (3,  3A,  and  3B).  The  nondimen- 

*  3|<  o  7 

sional  parameter  St  Pr  ’  is  plotted 
against  Re*j.  All  transport  properties 
were  introduced  at  the  reference  enthalpy 
and  static  pressure  for  each  station. 
Chemical  equilibrium  conditions  were 
used.  For  these  plots  duplicate  data 
points  represent  calculations  made  for 
each  of  the  last  two  common  time  reports 
for  each  run.  The  last  two  time  reports 
give  the  greatest  T  -  Tq  input  into  the 
parameter  <p  =  (T  -  TQ)/(Taw  -  Tq)  and 
therefore  have  the  least  error,  as  shown 
by  figure  7  (p.  10).  Figure  9  shows  that 
h  is  fairly  constant  with  time. 

In  each  part  of  figure  10  two  least- 
square  fits  for  the  data  are  plotted,  one 
where  the  constant  and  the  exponent  of 


Figure  11.  -  Exponent  and  correlation  constant  as  functions  of  .  ,  *  ,  ,  , .  ,  « 

area  ratio  for  two  correlation  equations.  where  the  constant  and  the  exponent  of 

the  Reynolds  number  were  both  allowed 
sk  skO  7  5kb 

to  vary  in  the  standard  equation  of  the  form  St  Pr  *  =  CRe^  ,  the  other  where  only 
the  constant  is  allowed  to  vary  and  the  exponent  b  is  held  at  a  fixed  value  of  -0.2.  The 
constants  and  the  exponents  are  shown  for  each  curve.  The  plots  also  have  the  line  rep¬ 


resented  by  the  equation  St '  Pr  '  '  =  0.  026  Re^_v"  ",  which  is  similar  to  the  simplified 

Bartz  equation  (ref.  12). 

In  order  to  give  a  better  statistical  interpretation  of  the  data  spread  around  the  much 
5k  *n  7  j|c  n  9 

used  St  Pr  =  CRe^  '  correlation,  the  standard  deviation  and  the  percent  deviation 
for  each  station  are  presented  in  the  preceding  table.  The  column  headed  1.  96  o  shows 
the  deviation  around  C  defining  the  95  percent  confidence  interval. 


*-0.  2 


,  which  is  similar  to  the  simplified 
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Over  the  Reynolds  number  range  covered  by  these  experiments  either  correlation 
(St*Pr*°‘  7  =  CRe*j"0' 2  or  St*Pr*®‘  7  =  CRe*jb)  seems  to  be  adequate  for  obtaining  heat- 
transfer  coefficients  at  a  given  station  for  the  hydrogen- oxygen  system  used. 

Figure  11(a)  is  a  plot  of  the  constant  against  area  ratio  for  an  exponent  of  -0.  2. 

These  data  indicate  that  no  one  value  of  C  will  correlate  the  data  for  all  axial  locations. 
In  fact,  the  constant  C  =  0.  026,  which  has  been  widely  used  to  predict  gas-side  heat- 
transfer  coefficients  for  nozzles,  overpredicts  the  coefficients  for  all  axial  locations 
except  the  chamber.  The  measured  heat-transfer  coefficients  at  the  throat  were  approx¬ 
imately  42  percent  less  than  those  predicted  by  the  equation  St*Pr*®‘  7  =  0.  026  Re*j  '  . 
Figure  11(b)  shows  plots  of  the  exponent  and  the  constant  against  area  ratio  where  the 
curve  fit  of  the  data  allowed  both  to  vary. 

From  hydrogen- oxygen  data  presented  in  reference  5  for  a  chamber  pressure  of 
150  pounds  per  square  inch  absolute  and  a  mixture  ratio  of  5.  0  (16.  7  percent  hydrogen), 
it  was  possible  to  compute  values  of  the  throat  heat-transfer  coefficient  h.  This  com¬ 
putation  gave  a  band  of  h  from  0.  0014  to  0.  0010  Btu  per  square  inch  per  second  per  °R, 
which  compares  with  a  value  of  0.  00136  measured  at  the  throat  of  the  nozzle  tested 
herein  under  the  same  operating  conditions.  Thus,  the  level  of  the  heat-transfer  coef¬ 
ficients  at  the  throat  shows  reasonable  agreement  for  the  hydrogen- oxygen  system  at  low 
chamber  pressures. 

Data  taken  from  reference  3  for  three  nozzle  configurations  were  recalculated  using 
a  reference  enthalpy  input  of  all  transport  properties  and  fitted  to  the  equation 
St*Pr*^'  7  =  CRe*j~^'2.  The  resulting  C  values  at  the  throat  were  0.  019,  0.017,  and 
0.  023,  respectively.  Data  taken  from  reference  4  were  recalculated  in  a  similar  man¬ 
ner,  and  a  value  of  C—  0.018  at  the  throat  resulted.  The  curves  of  C  against  area 
ratio  A/A*  also  showed  trends  similar  to  those  of  figure  11  (i.  e. ,  values  of  C  high  in 
the  chamber,  low  in  the  throat,  and  increasing  again  at  the  exit).  A  mention  of  this  work 
is  given  only  to  say  that,  while  the  variation  of  C  with  A/A*  that  was  found  is  con¬ 
sidered  correct  for  the  nozzle  geometry,  propellant  system,  and  injectors  used,  the 
values  should  not  be  used  universally  for  calculating  the  heat-transfer  rate  to  other 
nozzles. 


CONCLUDING  REMARKS 

If  a  designer  desires  to  use  the  correlations  presented  herein  for  a  system  similar 
to  the  rocket  of  this  experiment,  he  must  work  with  the  thermodynamic  and  transport 
properties  for  chemical  equilibrium  as  computed  herein  (refs.  9  and  10).  All  properties 
must  be  inserted  at  the  reference  enthalpy  and  static  pressure  of  the  axial  station  where 
the  heat-transfer  coefficient  is  desired.  The  effect  of  using  different  properties  can  be 
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great.  Using  the  new  transport  properties  (ref.  10)  reduced  the  correlation  constant  C 
by  approximately  30  percent  over  those  obtained  by  using  transport  properties  from 
reference  13.  Transport  properties  given  in  reference  10  are  considered  to  be  more 
accurate  than  those  previously  used. 

The  values  of  the  correlation  constant  C  reported  herein  are  for  a  given  nozzle 
geometry,  propellant  combination,  and  injector  configuration  and  are  not  to  be  con¬ 
sidered  universal  for  rocket  heat-transfer  calculations. 

SUMMARY  OF  RESULTS 

An  experimental  investigation  was  conducted  at  the  Lewis  Research  Center  to  de¬ 
termine  the  hot-gas-side  heat-transfer  rates  in  a  rocket  nozzle.  Transient  temperature 
measurements  were  made  at  five  axial  locations  in  a  copper  heat -sink  nozzle  with  ex¬ 
pansion  and  contraction  ratios  of  4.  64  with  gaseous  hydrogen  and  liquid  oxygen  as  pro¬ 
pellants.  The  nozzle  was  operated  over  a  range  of  chamber  pressures  from  150  to 
1000  pounds  per  square  inch  absolute  at  fuel  percentages  ranging  from  11  to  17,  with 
most  measurements  made  at  approximately  15  percent  hydrogen.  The  following  results 
were  obtained: 

1.  The  data  at  a  given  location  in  the  nozzle  are  adequately  correlated  by  the  equa¬ 
tion  St*Pr*0,  7  =  CRe^"0,  2,  where  St*,  Pr*,  and  Re*j  are  reference  Stanton,  Prandtl, 
and  Reynolds  numbers,  respectively,  and  C  is  a  constant. 

2.  The  constant  C  varies  with  axial  location  and  is  less  than  0.  026  for  all  locations 
except  the  chamber. 

3.  The  constant  C  at  the  throat  was  42  percent  less  than  the  widely  used  value  of 
C  =  0.026. 

Lewis  Research  Center, 

National  Aeronautics  and  Space  Administration, 

Cleveland,  Ohio,  March  16,  1965. 
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APPENDIX  -  SYMBOLS 


A 

b 

c 

c 

G 


H 

Haw 

h 

k 


L* 


<£ 

Nx 

p 

Pr 

Q 

q 


Re 


d 


S 


cross-sectional  area,  sq  in. 

exponent  in  equation 
St*Pr*°-  7  =  CRe*jb 

constant  in  equation 
St*Pr*°-  7  .  CRe*b 

specific  heat  of  materials, 
Btu/(lb)(°R) 

proportionality  constant; 


enthalpy 

adiabatic  wall  enthalpy, 

Hs  ♦  fiP  <Htol  -  Hs) 

heat-transfer  coefficient, 
Btu/(sec)(sq  in.  )(°R) 

thermal  coefficient  of  conduc¬ 
tivity  of  material, 

Btu/(sec)(in.  )(°R) 

characteristic  length  of  rocket 
engine,  in. ,  V^A^ 

length  of  heat  flow  path  in  rod,  in. 

nondimens ional  parameter,  hX/k 

pressure,  lb/sq  in.  abs 

Prandtl  number 

quantity  of  heat,  Btu/ cu  in. 

heat  flow  rate  per  unit  area, 
Btu/(sec)(sq  in. ) 

Reynolds  number  based  on  di¬ 
ameter 

entropy,  Btu/(lb)(°R) 


St 

Stanton  number 

T 

temperature,  °R 

T 

aw 

«Haw  ps> 

t 

time,  sec 

h 

time  1,  sec 

42 

time  2,  sec 

V 

volume,  cu  in. 

w 

total  weight  flow  rate,  lb/sec 

WQ 

u2 

oxygen  weight  flow  rate,  lb/sec 

X 

distance  normal  to  heat  surface, 

in. 

I 

nondimens  ional  parameter,  X/Sf 

P 

material  density,  lb/cu  in. 

(X 

standard  deviation 

T* 

nondimensional  parameter, 
h^t/kpc 

nondimensional  temperature, 

<T  -  T0>/<Taw  -  T0> 

3 

average  value  for  three  circum¬ 
ferential  measuring  stations 
at  throat  of  rocket  (station  3) 

Subscripts: 

A,B 

circumferential  locations  at 
throat  (station  3,  120°  apart) 

c 

chamber 

min 

minimum 

orf 

orifice 

s 

static 

th 

throat 

20 


tot 

total 

1  o  0 

5  ’  instrumentation  stations 

4,5 

w 

wall 

X 

distance 

Superscript: 

0 

zero  burning  time  or  start  of 
burning 

*  reference  condition 
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